Accurate description of optical precursors and their relation to 
weak-field coherent optical transients 
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We study theoretically the propagation of a step-modulated optical field as it passes through a 
dispersive dielectric made up of a dilute collection of oscillators characterized by a single narrow- 
band resonance. The propagated field is given in terms of an integral of a Fourier type, which cannot 
be evaluated even for simple models of the dispersive dielectric. The fact that the oscillators have 
a low number density (dilute medium) and have a narrow-band resonance allows us to simplify the 
integrand. In this case, the integral can be evaluated exactly, although it is not possible using this 
method to separate out the transient part of the propagated field known as optical precursors. We 
also use an asymptotic method (saddle-point method) to evaluate the integral. The contributions 
to the integral related to the saddle-points of the integrand give rise to the optical precursors. 
We obtain analytic expressions for the precursor fields and the domain over which the asymptotic 
method is valid. When combined to obtain the total transient field, we find that the agreement 
between the solutions obtained by the asymptotic and the exact methods is excellent. Our results 
demonstrate that precursors can persist for many nanoseconds and the chirp in the instantaneous 
frequency of the precursors can manifest itself in beats in the transmitted intensity. Our work 
strongly suggests that precursors have been observed in many previous experiments. 

PACS numbers: 42.25.Bs, 42.50.Gy, 42.50.Nn 



I. INTRODUCTION 

A fundamental problem in classical electromagnetism 
is the propagation of a disturbance or pulse through a 
dispersive optical material, which is characterized by a 
frequency-dependent complex refractive index n(u). The 
moment when the field first turns on (the pulse "front" ) 
is directly related to the flow of information propagating 
through the material The theoretical formulation 
is straightforward for the case when the input field is 
weak enough so that the dielectric responds linearly to 
the applied field. The calculation of the propagated field 
E(z,t), assumed to be an infinite plane wave traveling 
along the z-direction, involves the evaluation of a Fourier 
integral that is crucially dependent on n(u>) (see Sec. HIT). 
The exact evaluation of the integral is impossible even for 
simple causal models of the dispersive optical material. 

The first real theoretical headway on the problem was 
made nearly a century ago by Sommerfeld and Brillouin 
(SB), studying a step- modulated input field, which has 
zero initial amplitude and jumps instantaneously to a 
constant value Ao. Many aspects of the solution are sim- 
ilar to those for other pulse shapes. As summarized in 
a more recent collection of their earlier papers [l[ , Som- 
merfeld and Brillouin were able to show that the front of 
the step-modulated pulse always propagates at the speed 
of light in vacuum c and hence the flow of information 
is relativistically causal. They also found that, after the 
front and before the field eventually attains its steady- 
state value, there exist two transient wavepackets, now 
known as the Sommerfeld and Brillouin precursors. The 
wavepackets in the propagated field arise from contribu- 
tions to the Fourier integral that are localized near com- 
plex frequencies, known as the Sommerfeld or Brillouin 



"saddle-points," described in Sec. [II] The concept of pre- 
cursors exists only over the time in which such localized 
contributions to the integral occur. During this time, the 
sum of the Sommerfeld and Brillouin precursors provides 
the leading asymptotic approximation to the transient 
field. Over the years, various researchers have corrected 
errors in the SB calculations as well as extending the 
work to related problems, as discussed in great detail by 
Oughstun and Sherman (OS) @. It has been suggested 
that precursors can penetrate deeper into a material [31, 
which may be of use in underground communications [3( 
or imaging through biological tissue 

In this paper, we resolve a substantial controversy that 
exists on the observability of precursors in the optical 
part of the spectrum. The controversy was grounded in 
the belief that precursors are an ultrafast effect, where 
it is difficult to measure directly the field transients 
[S IE H 0, B H, 03 • Under incident frequency and mate- 
rial parameter conditions described below, we calculate 
expressions for the precursors and we find that their du- 
ration can be long (in the range of nanoseconds) and 
that their duration is controlled by the inverse of the res- 
onance width. We also find that the fraction of the tran- 
sient part of the field that is associated with precursors is 
controlled by the absorption. Furthermore, we reproduce 
Crisp's formula for the total field without making an 
assumption (see discussion and Appendix) that is crucial 
for Crisp's derivation. 

Our calculations of the precursors and of their lifetimes 
assume the following material parameter and frequency 
conditions. The half- width at half-maximum of the ma- 
terial resonance 5 is narrow, the carrier frequency of the 
field oj c is nearly equal to the material resonance fre- 
quency ujq, and the density of oscillators, characterized 
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by the plasma dispersion frequency u) p , is small enough so 
that certain limits are satisfied, as discussed later. These 
assumptions result in a simplification of the Fourier in- 
tegral for the propagating field that allows us to carry 
through the saddle-point calculation explicitly. Our ap- 
proximations differ from the approximate theory (and 
from the numerical evaluation) of OS toward the calcu- 
lation of the saddle-points. The OS theory is suitable 
for materials characterized by a broad resonance and a 
high density of oscillators, namely the situations when 
Lu, p and 6 are of the order of ujq. This restricted range of 
parameters was originally considered by SB [l[ and used, 
to a large extent, by most other researchers investigating 
precursor behavior. The OS approximations were not in- 
tended for dilute materials with narrow resonance and 
lead to unphysical predictions if applied to this case. For 
example, they yield an unphysically large amplitude of 
the transmitted field in weak-field coherent optical tran- 
sient experiments, recently performed by Jeong et al. @. 



The paper is arranged as follows. In Sec. [Til we formu- 
late the theory and, in particular, the Fourier integral for 
the propagated field. In Sec. IIII1 we simplify the Fourier 
integral based on considerations of the magnitudes of 
the material parameters and we apply the saddle-point 
method to the simplified integral to obtain highly accu- 
rate analytic expressions for the saddle-points and for 
the Sommcrfcld and Brillouin precursors. We find that 
the envelope of both precursors is non-oscillatory and 
that they display a frequency chirp. When these fields 
are added to yield the total transient propagated field, 
we find that the field envelope oscillates as a result of 
the chirp. In Sec. IIV1 we derive explicit asymptotic con- 
straints on the material parameters that guarantee the 
accuracy of our approximations used in the simplifica- 
tion of the Fourier integral. We make a direct, exact 
evaluation of this integral in Sec.|V] This calculation can 
only make predictions concerning the total transient field, 
but not the individual precursor fields. It gives identical 
predictions to the saddle-point theory under conditions 
when the latter is valid, as shown in Sec. IVII It agrees 
exactly with the result of Crisp obtained with the aid 
of the slowly varying amplitude approximation (SVAA), 
where we evaluate our theory in the limit where local field 
effects are negligible to make this comparison. Thus, we 
demonstrate unambiguously that the the weak-field co- 
herent optical transients resulting from the interaction 
of resonant radiation propagating through a dilute gas of 
atoms {e.g., the On pulse of Crisp [ll|) consist of optical 
precursors that can persist for many nanoseconds. We 
conclude that precursors have been observed in several 
experiments over the past few decades, discussed in Sec. 
IVIII A method for calculating higher-order corrections 
to further improve the accuracy of our results is given in 
the Appendix fSec. IVIII[) . 



II. THEORETICAL FORMULATION 

The propagated field E{z, £), assumed to be an infinite 
plane wave traveling along the z-direction, is expressed 
as the real part of a Fourier integral Q 



E(z,t) = -Re 



AqO(t) 
2tt 
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where the incident field is given by E(z — + ,t) — 
AoQ(t)sm(u c t). In Eq. (l}, 

T = t-z/c (2) 

is the retarded time, O(t) is the Heaviside function and 



, ' zn(uj) 
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where n(uj) is the refractive index at frequency u) and 
z is the depth of the measurement point. Equation (fT]) 
is an exact solution to Maxwell's equations, in integral 
form, for a step-modulated field propagating through a 
dispersive dielectric. 

In the original work of SB [l[ , and for much of the later 
work including that of OS Q, the dielectric is modeled as 
a collection of damped harmonic oscillators that fills the 
half space z > (known as a Lorentz dielectric). This 
model assumes that the density of oscillators (propor- 
tional to uip) is not too large. For high densities, the local 
field about an oscillator has a substantial contribution 
from its neighboring oscillators. Such local-field effects 
can be taken into account using the Lorentz-Lorenz for- 
mulation of the complex refractive index, which is given 

by Ema 
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The standard expression for the refractive index that 
does not take into account local-field effects can be ob- 
tained by dropping the last term in Eq. ([5]). The refrac- 
tive index depends on the frequency analytically except 
at the complex frequencies at which the square root has 
a branchpoint. This occurs when the quantity under the 
radical is cither infinite, in which case a — 0, with roots 



cj± = -i6 ± J - - <5 2 , 
or when it is zero, with a — ujp 1 , and roots 
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the counter-rotating term is 



FIG. 1: The complex oj-plane showing the original integra- 
tion path (dashed line, offset vertically for clarity) and the 
deformed contour of integration (solid line, offset vertically 
for clarity). Also shown is the singular point at u — ujo — iS 
and the simple pole at lo — ujq. 



For the material parameters of interest, the four roots 
lie in the lower complex half-plane, symmetrically about 
the imaginary axis, with imaginary part — iS. They are 
connected by two branchcuts, as shown in Fig. [TJ 

In order to facilitate the evaluation of the complex in- 
tegral for the propagated field, we deform the original 
contour of integration (i.e., the real axis, slightly moved 
to leave the pole under it), to a semicircle of infinite ra- 
dius in the lower complex half-plane that connects the 
± real points at infinity. The value of the integral over 
this semicircle is zero. As the deforming contour cuts 
through the two obstructing branchcuts, it leaves behind 
two clockwise-oriented contour loops, C+ encircling the 
right branchcut and C_ encircling the left one. As shown 
in Fig. [TJ the loops, except for their orientation, are cho- 
sen to be mirror images of each other with respect to 
the imaginary axis. Furthermore, is chosen to pass 
through the Sommerfeld and Brillouin saddle-points (i.e., 
stationary points of the exponential in the Fourier inte- 
gral, see below) in the right half-plane. This implies that 
C_ also passes through the two saddle-points in the left 
half-plane. Within the regime of validity of the saddle- 
point method, the main contributions to the integrals 
are localized near these saddle-points. The pole residue 
contribution must be added to the field if the pole uj c is 
positioned outside C+. 

The integral over C_, representing the so-called 
counter-rotating contribution to the propagated field, can 
be efficiently represented as an integral over C+ through 
a change of the variable of integration u — > — ZU (bar indi- 
cates complex conjugate) and using the symmetry of the 
exponent, 



%l>(-u) = ip(ui). 



We thus obtain 



-du> = — 



c u + w c 



-dm. 



Inserting Eq. @ in Eq. ([1}, we obtain 

E(z,t) = E rot (z,t) + E crot (z,t) + £, C+ E C , 
where the rotating term is given by 

, . / A e(r) f , N 

E rot (z,t) = -Re [ °. 1 ; / du 



2tt j c u. - 



(8) 
(9) 
(10) 

(11) 



E crot (z,t) = Re I — / du 



the pole contribution is 



2TT J(j ui + ui c 



E c = A e(T)e' aaZ / 2 sin ( w r 



oj()An r 



(12) 



(13) 



which represents the main field, and £c + = or 1 accord- 
ing to whether C + encloses the pole or not (if the pole is 
enclosed, the main field is contained in E rot ). The sum of 
the rotating term and the third term in Eq. (|10[) remains 
constant under deformations of the contour C+ . 

The expression for the main field (|13|) displays expo- 
nential attenuation as a function of propagation distance, 
which is governed by the absorption coefficient ao at fre- 
quency u) c = u/q. It is defined through the relation 



a 



2aj n l (u ) 



where 



rii(ui ) = Im[n(w )]- 



(14) 



(15) 



The value of the refractive index at the resonant fre- 
quency is given by, 



n(cj ) = 



l + im/3 
1 — im/6 ' 



or, in polar form, 



n(u) ) = 2 



m 2 + 9 iu^-i 



m? + 36 



where 



i[tan _1 9m/(18-m 2 )]/2 
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Thus, 



ni(uJo) = Im[n(u> )} 



(16) 
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and 



n r (u Q ) 



Re{n(u )] 
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(20) 
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The main signal also experiences a ^-dependent phase 
shift arising from the real part of the refractive index, 
where An r — n r (uj ) — 1. As discussed in Sec. I VIII our 
derivation of expressions for optical precursors allows for 
large absorption coefficients and Eqs. (fT^|) and ([2H)l must 
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be used. On the other hand, a simplified expression for 
the refractive index can be obtained, and a connection 
to other treatments of optical pulse propagation can be 
made when the absorption is small (ao <C ujq/c). In this 
case, 



i( w o) 



n r (uj ) ~ 1, 



a 



2c~5' 



An r ~ 0. 



(21) 
(22) 

(23) 
(24) 



Here, we see immediately that ao scales with w^ and 
inversely with S. 

Using Eqs. (TTT jl -lfTB" ! in the expression for the field (fT0| . 
we calculate the precursors from contributions from only 
the Sommerfeld and Brillouin saddle-points in the right 
half-plane; they include the contributions from their sym- 
metric counterparts in the left half-plane through the sec- 
ond integral. Henceforth, references to saddle points or 
a branchcut are to the ones in the right half-plane. 

In order to perform the calculation of the saddle-points 
explicitly, and thus obtain an explicit expression of the 
precursor fields, we focus on an asymptotically large ma- 
terial resonance frequency u)q. In this limit, the material 
is dilute (lo p <C wo) and narrowbanded (5 <C Wo). In the 
scale of wo, the branchpoints © and ([7]) of the right half- 
plane (w + and lu' + , respectively) collapse asymptotically 
to the "singular point" loq — iS, as illustrated in Fig. [1] 
The more precise asymptotic formula for the midpoint of 
the collapsing branchcut (not needed in our calculation) 
is uj — i<5+(tjp/12 — 6 2 /2)/ojq. The length I of the branch- 
cut, i.e., the difference lo' + — w + , is given asymptotically 
by 



2w 



< w . 



(25) 



We make the singular point the center of a new fre- 
quency variable, denoted by w* and defined by 



UJ 



ujq — iS + w* 



(26) 



We seek parameters (frequency, material, depth, retarded 
time) for which the values of w» at the Sommerfeld and 
Brillouin saddle-points are at a scale that is intermediate 
between the large material frequency loq and the small 
length I of the branchcut. We require 



2ujr 



•C w* <C wo- 



(27) 



As a result, the saddle-points view the branchcut as a 
point. Furthermore, for near resonance excitation (lj c ~ 



wo), the saddle points are disproportionately farther away 
from the center frequency — loq+iS of the counter-rotating 
term than from the center frequency wo + iS of the rotat- 
ing terms. Using these facts allows us to obtain explicit 
expressions for the precursors and for the transient field, 
as discussed below. 



III. CALCULATION OF THE PRECURSORS 

In order to calculate the integrals that give the field, 
we seek to isolate the dominant terms in the exponent 
tp and, in particular, in the expression for the refractive 
index. The value of a in Eq. expressed in terms of 
the new frequency variable w», is 



2w w* + luI + -Wp. 



Following our scaling, 



2w>ow* 



and the term ojp/a in Eq. Q satisfies 



v 
a 



< 1. 



Defining 



4w w* ' 

we write the refractive index as 

n = 1 — m + d 

where the error term d satisfies d 
Eq. (J32J) in Eq. © for ip to obtain 



(28) 



(29) 



(30) 



(31) 



(32) 

0{ri[). We insert 



-iujr — i—uj{ni — d). 



(33) 



where we have used the retarded time (Eq. ©). 

We insert the change of variable (|26|) into this expres- 
sion and arrange the terms into three groups, as indicated 
here with the aid of braces 



%j) — {-iuj T — St} + j-zw*T — i-wonxj (34) 

f z z z ~\ 

+ < i—ujQd 5(ni — d) — i— w*(ni — d)> . 

I c c c ) 

The two terms of the first group are independent of the 
variable of integration and give the leading amplitude 
and phase contributions to the integral. Our theory ap- 
plies to parameters (to be identified below) for which the 
second group is dominant, allowing for the third group, 
labeled iprem (mnemonic remainder) to be neglected in 
the calculation of the saddle-points. The second group is 
further simplified by introducing the rescaled frequency 
variable r\, defined by 



w* 
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ZV, 



(35) 
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where q is given by 



(36) 



Using these notations, the phase is given by 

ij) = {-iuoT - 5t} + {-iq(rj + ii^ 1 )} + ip rem . (37) 

Inserting this expression into Eqs. (Ilip and (fl"2")) and 
changing the variable of integration to rj, we obtain 



E rot (z,t) 



E C rot(z,t) = Re 



-Re 



Ai 
2ir 



Ai 
2ir 



c+(n) 



T] — IO 



-dr\ 



c+(v) 



1] — ia + 2(w c r / 1 q) 
Ax = A e{T)e- 5T -* TUJO , 



(38) 

ch] 
(39) 
(40) 

(41) 



where C+ (rj) is the image of the contour C+ in the rj plane 
as shown in Fig. We note that formulae (|3"8"|) and (|3U)) 
are exact, in spite of their derivation being guided by 
asymptotic considerations. 



C + (ri) 
B(il=-1) 



S(r,=1) 




Re r\ 



FIG. 2: The complex ry-plane showing contour of integration 
Cd (dashed line) and the short stretches along d that give 
the dominant contribution to the integral. The contour Cd is 
a deformation of the contour C3 (not shown for clarity) . The 
solid line is a circle of radius 1 centered at the origin. The 
shaded regions are the portion of the 77-plane for which the 
real part of the exponent of the integrand is negative. 

We now turn to our calculation of the precursors, which 
implements the following approximations. 

1. We neglect the terms of iprem altogether, i.e., we 
set -0rem = (see discussion in section HV)) . The 
saddle-points are then approximated as the station- 
ary points of 77+77^ 1 , namely, 77 = 1 for the Sommer- 
feld precursor and r/ = — 1 for the Brillouin precur- 
sor (see Fig. [5]). The saddle-point approximation 
of the integrals is straightforward and given below. 



The omitted terms give only higher order contribu- 
tions to the phase and amplitude of each precursor. 
A method for obtaining higher-order corrections is 
discussed in the Appendix. 

2. The counter-rotating term of the field E crot is ne- 
glected. In our case of near-resonant excitation and 
with our scaling, the relative error introduced is of 
order q/(TLOo) = (lOp/ujq)^ zj (cr) <C 1 (see relation 
(j60|) below), following the fact that = q/r at 
the saddle-points. 

As a result of these approximations, the propagated field 
is given as 



E(z,t) = -Re 



M 

2ir 



c+(v) ■n - 



-dr) +e(<7-l), 
(42) 

where C+(rj) is the unit circle in the rj plane (see Fig. [2]). 

The saddle-point method for evaluating the integral in 
Eq. (|42[) requires large values of q, the relative error of 
its approximation to the value of the integral being of 
order 1/q. In order to apply the method, we deform the 
contour of integration C + (ri) to the contour Cd (see Fig. 
(2), oriented clockwise, that passes through the saddle- 
points rj = ±1 cutting the real axis at angle tt/4. The 
value that the real part of the exponent assumes at the 
saddle-points is maximal, in comparison to its values in 
the shaded region shown in the figure. The contour of in- 
tegration Cd passes through the saddle-points and stays 
in the shaded region. In the limit of large q, exponen- 
tial decay in the shaded region makes the contribution of 
the part of the contour close to the saddle points domi- 
nant. The main contributions to the integral arise from 
the two short stretches of the contour Cd in the neigh- 
borhood of the saddle points, which have been chosen in 
the steepest descent direction (-7r/4 angles with the real 
axis) where the exponential is purely real and the inte- 
grands are approximated by Gaussians. Large values of 
q localize the Gaussians at the saddle-points. Thus, the 
length of the stretches tends to zero as q increases. The 
contributions from the two saddle-points to the rotating 
term of the field, obtained from the exact calculation of 
the Gaussian integrals, are 



E s (z,t) = Re 



E B (z,t) = Re 



^ ie -*(2?-f) 
_2(nq)^ (1-ia) 

2 {nqY (1 + ia) 



(43) 



(44) 



with the subscript S (B) for the Sommerfeld (Brillouin) 
precursor field. Inserting Eq. ([50]) into these obtains the 
two precursor fields 



E s (z,t) = Re 
E B (z,t) = Rc 



iA 6(r)e 



-8r-t(2q-f) 



2(-Kq) 2 (1 - ia) 

iA Q(T)e- ST + i( .^-^ 

1 1 

2(nq) 2 (l + ia) 



(45) 
(46) 
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The precursors display rapid oscillations at a frequency 
close to lu 0i modulated by a complex- valued envelope. 
Both precursor envelopes have the same modulus 



A e(r)e 



-St 



2(Trq{l + a 2 )y 



(47) 



The precursors decay exponentially with time constant 
1 /S, supporting our statement that they persist for a time 
determined by the resonance half-width, which can be in 
the nanosecond time scale for a dilute gas of cold atoms 
@, for example. The precise value of the frequency of 
the Sommerfeld and Brillouin precursors is determined 
by taking the derivative of the phase TajQ+2q with respect 
to r in Eqs. (|45|) and (|46j) . respectively, yielding 



T 



q 

UJQ . 

T 



(48) 



Note that the precursors frequencies are equal to the real 
part of the respective saddle points in the complex em- 
plane. Figure [3] shows the envelope and frequencies of 
the precursors using the materials parameters of the ex- 
periment of Jeong et al. [i| with ljq = 2.5 x 10 15 s -1 , 
Ljp = 3 x 10 9 s _1 , 5 = 3 x 10 7 s _1 , lj c = uj , but with 
a longer medium length z=20 cm. It is seen that the 
envelop persists for many nanoseconds and that the pre- 
cursor frequencies are within a few hundred MHz of the 
resonance within a few nanoseconds. 

The total transient field is approximated by the sum 
of the two precursor fields 



E T (z,t) =E s (z,t) + E B (z,t). 



(49) 



Because the individual precursors are chirped, Et will 
display oscillations at the beat frequency los—ujb = 
whose period increases with r (the beat frequency de- 
creases with t) . The envelope decays with time constant 
1/(5 as do the individual precursors. The expression for 
the envelope of the total precursor field is given by 



A T {z, t) = 2A s ,b{z, t)cos(2g - - 



tan V) 



(50) 



where q and a are given by (|3G[) and (|4"Tj) . respectively. 
The pole, located on the imaginary axis at ia, starts at 
the origin when r — and moves up the imaginary axis as 
time increases. Its contribution, i.e., the main field E c , is 
negligible compared to the precursor field when the pole 
is in the shaded region. Outside the shaded region, the 
main field is dominant. 



IV. PARAMETER ANALYSIS 

The determination of the range of parameters for which 
our calculation of the precursors is accurate follows from 
assumptions we have made, which we now summarize. 





1 .0 ■ 






0.8 ■ 




o 

< 








0.6 ■ 




N 






CO 

< 


0.4 ■ 
0.2 ■ 






0.0 ■ 





N 

X 



CNJ 



3^ 

I 

CO 

3°' 



300 
200 
100 


-100 
-200 
-300 




FIG. 3: a) The envelope and b) frequency of the Sommerfeld 
and Brillouin precursors field with the medium parameters of 
Ref. H and z = 20 cm. 



1. The requirement of small n\ in the expression for 
the refractive index (|3"2"|) gives 



«1, 



2lu 



(51) 



(52) 



This constraint is identical to the left relation (|27p . 
Inserting lu* — q/r, we obtain, after simple algebra, 



lj ti / z 
— < \ — ■ 

LOq V CT 



(53) 



2. The requirement of large q yields the constraint 

c 



LO p > 2 



(54) 



which can be obtained from the definition of q. 

3. The requirement of the asymptotic vanishing of the 
second term of the ip rem (third group of terms of 
Eq. O), 



: <5(rii -d) < 1, 



(55) 
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which, thus, does not contribute to the leading or- 
der of the integral, partly justifying our neglecting 
of iprem- The requirement, expressed in terms of 
material parameters, is 



for ni <C 1 and dominance of the second group in ip; and 



ujpS 
w 



ZT 
— < 1. 

C 



(56) 



4. The dominance requirement constrains the remain- 
ing (first and third) terms of ip rem to be signifi- 
cantly smaller than the terms of the second group 
along the contour of integration. Since the two 
terms of the dominant group have comparable mag- 
nitudes, it suffices to make the comparison with 
only the second term in the dominant group. The 
ratios of the first and third terms of iprem by the 
second term of the dominant group are, respec- 
tively, 



w.: 



2wow* 

W* 

w 



Up 

2w 



UJp / CT 
OJ \J Z 
UJn I Z 

ii V CT 



«1, 
< 1. 



(57) 
(58) 



These relations are identical to the scaling (left and 
right, respectively) of relation (127| . which is thus 
satisfied automatically as a result of the require- 
ment. 

By neglecting these two terms of ip rem , while falling 
short of requiring their asymptotic vanishing, we 
lose a phase term in the leading order expression 
for the precursors. That the error made is only in 
the phase follows from the fact that, after setting 
the second term of iprem equal to zero, the saddle- 
points are real (in the rj variable) and the exponent 
is purely imaginary. The dominance requirement 
guarantees that the error in the phase is of higher 
order compared to the phase correction from the 
second group in %p that produces the chirp. In terms 
of physical insight gained by our result, tolerating 
this error is preferable to further constraining the 
material parameters. The dominance requirement 
also guarantees there is no other leading order error 
in the application of the saddle-point method. 

Collecting the independent constraints leaves us with 



and 



W 



UJr,5 I ZT 

— \ — < 1. 

w V c 



(61) 



(62) 



for the asymptotic vanishing of the second term of iprem- 
When these constraints are satisfied, (a) the precur- 
sors exist at the specified retarded time {i.e., the calcu- 
lation of the transient field as the sum of saddle-point 
contributions applies) and (b) our calculation of the pre- 
cursors is accurate. When some constraint is violated, 
one of these statements may not be true. Assuming suit- 
able fixed frequency, depth and material parameters, con- 
straints (159")) and (|6ip can be satisfied only past a (usually 
short) retarded time r. The constraints are satisfied in a 
time-range beyond this, until, for time large enough, con- 
straint (|60)) is necessarily violated and our method loses 
accuracy. 

The upper time- limit of validity of our constraints may 
be overshadowed by the additional practical constraint 
that St must be fairly small for the precursor to be ob- 
servable. When 



s*_ 

,2 



z 

CT 



(63) 



the real exponentials multiplying the amplitudes of the 
precursor and of the main field, respectively, are equal. 
To be solidly in the regime where the precursor dominate 
over the main field, we require 



S 2 z 
-» < — 



(64) 



Generally, when relation (|59j) is comfortably satisfied, 
but some other constraint fails, we expect that the pre- 
cursors exist, but our calculation loses accuracy. To gain 
accuracy, we apply a corrective scheme described in the 
Appendix. 



EXACT EVALUATION OF THE INTEGRAL 
FOR THE APPROXIMATE FIELD 



2 J — < w E 



(59) 



for large q (the requirement of a large value of q is mod- 
est - the saddle-point method already gives a quite good 
approximation to the value of the integral for a value of 
q of 3 or 4); 



w 



(60) 



Equation (|42j) gives the approximate propagated field, 
in which ip rem and the counter-rotating terms have been 
ignored (the same approximations used to obtain the 
expression for the precursor fields). In order to make 
an exact calculation of the integral in this equation, we 
change to an angle variable of integration defined through 
rj = cxp[i(p + 7r/2)]. We thus obtain 



E{z,t) = Re 



iAi 
~2T 



1 - a&- % c 



dp 



-Q(a-1)E C . (65) 
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When the integrand of (|65)l is expanded in a series 
of powers of a < 1, each of the integrals in the series 
represents a Bessel function and the total field is given 
by 



E(z,t) = Re 



iAiJ2a k J k (2q) 



k=0 



(66) 



where Jk is the Bessel function of order fc. The ap- 
proximate transient field is obtained when the pole- 
contribution to the field (fT3|) , representing the main field, 
is subtracted from Eq. (f6"6) . While the ensuing expres- 
sion for the transient field is exact in the limit considered 
here, it does not allow separating the two precursor fields 
and does not even make a statement about the existence 
of individual precursors. This is due to the fact that the 
definition of the precursors is tied to the application of 
the saddle-point method in the calculation of the field. 

A similar calculation is possible when a > 1. In this 
case, the expansion is in powers of er -1 and the final 
formula is again a series of Bessel functions. Although we 
only have treated the case of a resonant field (u> c = too), 
the result for a near resonant field is obtained by inserting 
an imaginary part with the frequency difference in a and 
expanding in powers of a or ct _1 according to whether 
cr| is less than or greater than unity. 

In the parameter range of negligible local field effects 
(ao *C ojq/c) , an alternative formula for a is 



a z/2' 



(67) 



Interestingly, when this expression for a is used, Eq. (|66|) 
is identical to that found by Crisp who used the slowly 
varying amplitude approximation (SVAA) (see Eq. (29) 
of Ref. [HI). We derive Eq. (|66[) only with assumptions 
about the material properties; we make no assumption 
about the slowness of the variation of the electromag- 
netic field. Our formula (]66p is still valid, even when the 
relation ao <C loq/c is violated, as long as a is defined by 
Eq. gl|). 



VI. RESULTS 

For the first time, we have derived analytic expres- 
sions for optical precursors for a material with a nar- 
row resonance and a low oscillator number density (small 
plasma frequency) . Our formulae are valid within explic- 
itly specified sub-ranges of the space-time range of the 
existence of precursors. The latter consists of the range 
of points in space-time over which, field contributions to 
the main Fourier integral (Eq. |T|)) are localized at iso- 
lated saddle-points of the complex frequency plane. Our 
precursor theory limits itself to such saddle points that 
are sufficiently close to the resonant frequency for the 
counter-rotating field contributions to be negligible and 



sufficiently far from it to allow approximating the branch- 
cut in the refractive index by a singular point. These 
restrictions apply to the long tail of the precursors, i.e. 
after the narrow front of the transient wave has passed. 
The saddle-points involved are isolated. The front (it dis- 
plays degenerate or near-degenerate saddle-points) and 
the very early tail, both requiring the counter-rotating 
contributions, are not addressed in this study. The exact 
expression of the field is written in the form of Eqs. (|38[) . 
(fMjl . pi)) . (JUJ). The approximate field is given by Eq. 
(jl2"[) . on which the saddle-point method is performed to 
yield the precursors (jl5|) and (pr6"|) . The space-time and 
material parameter constraints that define the range of 
validity of the derivation of the approximate field are 
given by the relations (SB, dMJ) and (HU). The ex- 

act evaluation of the approximate field ([12")) , in terms of 
Bessel functions is given by Eq. ([66]). 

Our method can handle comparatively high number of 
oscillator densities, for which the condition ao <C loq/c 
and, hence, the a priori assumption of the SVAA (slowly 
varying amplitude approximation) are no longer valid. 
Clearly, local field effects play a significant role in the ex- 
pression of the main field. On the other hand, local field 
effects are still negligible in the expressions for the precur- 
sors. Indeed, the separation between the saddle-points 
and the resonant frequency is sufficiently high and the 
value of the refractive index at the saddle-points remains 
essentially unaffected. One verifies a posteriori, that the 
SVAA holds in the derived formulae for the precursor 
fields, both in time and in space. Indeed, the separation 
of the scales of the carrier and beat frequencies is guaran- 
teed by the relation dq/dr <C loq (equivalent to constraint 
([61 [) ) and slow amplitude variation in time is guaranteed 
by d <C luq . The small variation of amplitude in space is 
guaranteed by relation (1/ (2q))dq/dz <C 2dq/dz. (equiv- 
alent to 1 <C <?). The left side of this relation is exponen- 
tial attenuation, obtained from bringing the rational at- 
tenuation to the exponent as a logarithm and taking the 
spatial derivative. Our results are consistent, in the sense 
that the sum of precursors agrees with the transient field 
obtained from the exact evaluation of the approximate 
field in Eq. (H2) (see Fig. [|]). Our exact expression, in 
terms of a series of Bessel functions, is valid in the higher 
density regime as well and agrees with Crisp's formula 
when restricted to low densities (see discussion below). 

We now give an example of the predictions of the pre- 
cursor theory and compare the results to the exact cal- 
culation of the field for the dilute narrowband dielectric. 
In particular, we use the parameters of the experiment 
of Jeong et al. 9] with lu = 2.5 x 10 15 s"> p = 3x 10 9 



3 x 10 7 s~ 



coq. We first consider the case 



when the medium length is 0.2 cm, the value used in the 
experiment. For these parameters, relation ([59[) guaran- 
teeing a large value of q, becomes an equality (q takes the 
value q — 1) at r rs 67 ns, indicating that it takes on the 
order of hundreds of nanoseconds to satisfy the require- 
ment. Figure 2^ compares the transient field envelopes 
for the two theories. While there is some discrepancy at 
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FIG. 4: The envelope for the total transient field with the 
medium parameters of Ref. [9j] and a) z = 0.2 cm and b) 
z = 20 cm. The solid line shows the predictions of the exact 
theory and the dots show the predictions of the asymptotic 
(precursor) theory. 



shorter times, the error is less than 25% for times greater 
than 30 ns, and very small at 67 ns, indicating that con- 
dition (|59")) is rather conservative in this case. Thus, the 
primary contribution to the transient field is from the 
saddle points and hence it is reasonable to conclude that 
the experiments of Jeong et al. observed optical pre- 
cursors. We note that they found that the exact theory 
agrees very well with the experimental observations. 

Figure UJd compares the two theories for a medium 
length 100 times longer than that used in the experiment 
(z—20 cm), where it is seen that the agreement between 
them is excellent. The increasing-period oscillations in 
the transient field is clearly evident and consistent with 
our discussion above. In this case, q = 1 at r ~ 6.7 ns; 
the good agreement at as low a value as q = 1 again 
points to the conservatism of the condition (|59|) . 



VII. DISCUSSION 

Crisp found that the step-modulated input field evolves 
toward a so-called 07r pulse whose envelope oscillates. 
(Note that a step-modulated incident field inherently 
violates the SVAA, yet our solutions are valid for this 



situation.) He showed that the pulse area of the total 
field approaches zero, which is known as a 07t pulse in 
the quantum optics community. Such pulses have been 
studied experimentally by a number of groups, beginning 
with the observation of Rothenberg et al. [14j , later work 
demonstrating 07r-pulse 'stacking' [TEj . and more recent 
work 0, H3, Crisp posed the question of whether 
these weak-field coherent optical transients (the 07r pulse) 
are a manifestation of optical precursors and answered it 
in the negative without mathematical proof. Certainly, 
Crisp explained these oscillations as an interference be- 
tween the parts of the pulse spectrum above and below 
the atomic resonance frequency, where the central part 
of the spectrum is eaten away as the field propagates 
through the material. He did not associate these fre- 
quency components with the frequencies of the Sommer- 
feld and Brillouin precursors, reasoning that precursors 
are an ultrafast effect, which would violate the assump- 
tion of a slowly-varying amplitude, and thus must be 
precluded from the SVAA formalism. Our precise math- 
ematical analysis proves conclusively that Crisp's conclu- 
sion is incorrect. It follows clearly from our analysis that 
the oscillations in the envelope of the 07r pulse is the re- 
sult of the interference of the Sommerfeld and Brillouin 
precursors. 

Avenel et al. [1] , on the other hand, first suggested that 
the coherent optical transients predicted by Crisp and ob- 
served by Rothenberg et al. [14j are a manifestation of 
optical precursors and that the time scale for the precur- 
sors can be very long (of the order of nanoseconds) for 
a material with a narrow resonance. However, they did 
not provide a mathematical justification for their claim. 
In later work, Varoquaux et al. [19j attempted to make 
their claim precise by solving Eq. |T]) using an asymp- 
totic method. They found a solution to the integral only 
for frequencies well above the frequency of the material 
resonance luq (see the discussion near the end of their 
Sec. IV. D.). They predict that the envelope of the Som- 
merfeld precursor contains oscillations similar in form to 
that predicted by Crisp, which is the same as our exact 
solution (Eq. ([66]) ). They were not able to identify a Bril- 
louin precursor. It is not surprising that they failed to 
obtain an accurate prediction concerning the precursors 
because the saddle points are located close to o>o (with 
respect to the scale of ujq), yet their approximate solu- 
tion to the integral only accounted for the contributions 
to the integral at much higher frequencies. Our calcu- 
lation identifies all the saddle-point contributions to the 
integral and places the Avenel et al. conjecture on a firm 
theoretical foundation. 

In addition to the calculation of the precursors, we 
calculated the total propagated field through the exact 
evaluation of the simplified Fourier integral (|4"2"|) . The 
calculation gives the total propagated field as a series of 
Bessel functions. This result was first obtained by Crisp 
[ll| under the additional assumption of the slowly vary- 
ing amplitude approximation (SVAA). In the SVAA ap- 
proach, the wave equation is first simplified by assuming 
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a slowly- varying amplitude, then approximations about 
the material are invoked, and a solution is thus obtained. 
The use of a step-modulated nature of the initial field in 
the context of the SVAA has raised questions. The exact 
agreement of Crisp's formula with our results (evaluated 
for a low density of oscillators), which makes no assump- 
tion of slowly varying amplitudes, demonstrates that, in 
spite of the initial discontinuity, the slowly-varying as- 
sumption is superficial for a weak-field and a narrow- 
resonance dilute medium. In fact, it can be shown that 
the solution to the SVAA equations for these conditions 
also is a solution to the full wave equation. Our deriva- 
tion shows the correct way to extend Crisp's formula to 
the high-density regime. 

Finally, we remark that our work also has implica- 
tions for very weak 'quantum' fields. Very recently, Du et 
al. [13, [ll| have shown that precursors can be observed 
on long time scales in correlated bi-photon states. A fol- 
low up study considering the propagation of a classical 
field through a similar medium has also been presented 
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where r\ = 1 + 771 , and 70 and 71 are the values of ip rem 
and its first derivative at the Sommerfeld saddle-point, 
respectively. After letting rj~ 1 = (1 + 771)^ 1 ps 1— rji+r/f, 
we obtain for the exponent in 



(70) 



( , 1 1 ~ fr, , 2 , 70 . 71 

- q(v + -) « -9 2 + r)x + — 1 — m, , 

V \ 1 1 

for the stationary point 

7i 
2q 

for the exponent maximum 



2<Z- 7 o + 



(71) 



(72) 



and the second derivative at the stationary point 

- 2q. (73) 

In order to insert the corrections to the precursor field 
(|45| . we 

1. Adjust the amplitude by multiplying the field by 
the factor 

„ — ro+7i/ 4 9 



VIII. APPENDIX: HIGHER-ORDER 
CORRECTIONS 

When the relation 

Tprern < q{V + T]' 1 ) (68) 

does not apply comfortably, we cannot neglect iprem al- 
together in determining the saddle-points. Instead, we 
correct our calculation of each precursor by retaining 
the linear Taylor approximation of iprem about the corre- 
sponding saddle point. At the Sommerfeld saddle-point 
rj = 1, for example, the linear Taylor approximation of 
4> rem is given by 

Iprem = 70 + JlVl, (69) 



which inserts the correction of the exponent max- 
imum (the second derivative of the exponent at 
the stationary point whose square root enters the 
fromula for the saddle-point contribution remains 
unchnaged in corrected exponent). 

2. Perform the replacement 

in the denominator. 

The calculation is similarly straightforward for the Bril- 
louin precursor. The corrective procedure maybe iterated 
by using the updated saddle-point as the base point. 
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